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Abstract. Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., gamma- 
ray bursts (GRBs), active galactic nuclei (AGNs), and microquasars commonly exhibit power-law emission spectra. Recent 
PIC simulations of relativistic electron-ion (or electron-positron) jets injected into a stationary medium show that particle 
acceleration occurs within the downstream jet. In collisionless, relativistic shocks, particle (electron, positron, and ion) 
acceleration is due to plasma waves and their associated instabilities (e.g., the Weibel (filamentation) instability) created 
in the shock region. The simulations show that the Weibel instability is responsible for generating and amplifying highly non- 
uniform, small-scale magnetic fields. These fields contribute to the electron's transverse deflection behind the jet head. The 
resulting "jitter" radiation from deflected electrons has different properties compared to synchrotron radiation, which assumes 
a uniform magnetic field. Jitter radiation may be important for understanding the complex time evolution and/or spectra in 
gamma-ray bursts, relativistic jets in general, and supernova remnants. 
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INTRODUCTION 

Shocks are believed to be responsible for prompt emis- 
sion from gamma-ray bursts (GRBs) and their after- 
glows, for variable emission from blazars, and for par- 
ticle acceleration processes in jets from active galactic 
nuclei (AGN) and supernova remnants (SNRs). The pre- 
dominant contribution to the observed emission spectra 
is often assumed to be synchrotron- and inverse Comp- 
ton radiation from these accelerated particles ifTsi 
It is assumed that turbulent magnetic fields in the shock 
region lead to Fermi acceleration, producing higher en- 
ergy particles JUITt]. To make progress in understanding 
emission from these object classes, it is essential to place 
modeling efforts on a firm physical basis. This requires 
studies of the microphysics of the shock process in a self- 
consistent manner lIlHfTgll . 

New Numerical Method for Calculating 
Emission 

The retarded electric field from a charged particle 
moving with instantaneous velocity p under acceleration 
j3 is obtained . After some calculation and 



simplifying assumptions the total energy W radiated per 
unit solid angle per unit frequency can be expressed as 
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Here, n = R(t ) / |R(t' ) | is a unit vector that points from 
the particle's retarded position towards the observer. The 
first term on the right hand side, containing the veloc- 
ity field, is the Coulomb field from a charge moving 
without influence from external forces in eq. 2.4 (0]. 
The second term is a correction term that arises when 
the charge is subject to acceleration. Since the velocity- 
dependent field falls off in distance as R^^, while the 
acceleration-dependent field scales as R^^, the latter be- 
comes dominant when observing the charge at large dis- 
tances {R^ 1). The choice of unit vector n along the di- 
rection of propagation of the jet (hereafter taken to be the 
X-axis) corresponds to head-on emission. For any other 
choice of n (e.g., = 1 //), off-axis emission is seen by 
the observer The observer's viewing angle is set by the 



choice of n (n^ + n^ 



1). 



TABLE 1. Seven cases of radiation 
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RADIATION FROM TWO ELECTRONS 

In the previous section we discussed how to obtain the 
retarded electric field from relativistically moving parti- 
cles (electrons) observed at large distance. Using eq. 1 
we calculated the time evolution of the retarded electric 
field and the spectrum from a gyrating electron in a uni- 



form magnetic field to verify the technique used in this 
calculation. We have calculated the radiation from two 
electrons gyrating in the x — y plane in the uniform mag- 
netic field Bz with Lorentz factors (7= 15.8,40.8) (Case 
P in Table 1) lfT2l[T3l [T?l. We have very good agreement 
between the spectrum obtained from the simulation and 
the theoretical synchrotron spectrum expectation from 




FIGURE 1. Summary of six cases with jet velocity (Cases A - F) (for Case P, see Nishikawa et al. 814ll ) 



eq. l.iof^. 

In order to calculate more realistic radiation from rel- 
ativistic jets we included a parallel magnetic field (fix)- 
Relativistic jets are propagating along the x direction. 
Table 1 shows six cases including the previous case P 
(first row) |14]. The jet velocity is 0.99c (except Case 
B). Two different magnetic field strengths are used. Two 
electrons are injected with two different perpendicular 
velocities (Cases A - F). The maximum Lorenz fac- 
tors, 7max = {(1 - (Vj2 + VIj)/^^}"^''^ are calculated for 
larger perpendicular velocity. The critical angles for the 
off-axis radiation is calculated with 0r = ' . 

Figure 1 shows the summary of the six cases. Tra- 
jectories of the two electrons are shown in the left col- 
umn (red: larger perpendicular velocity, blue: smaller 
perpendicular velocity). The two electrons propagate 
from left to right with gyration in the y — z plane (not 
shown). The gyroradius is about 0.44A (A = 1: the sim- 
ulation grid length) for the electron with a larger per- 
pendicular velocity (Case A). The radiation electric field 
from the two electrons is shown in the middle col- 
umn. The spectra were calculated at the point {x,y,z) = 
(64,000,000.0,43.0,43.0) shown in he right column. 
The seven curves show the spectrum at the viewing an- 
gles 0° (red), 1 ° (orange), 2° (yellow), 3° (moss green), 
4° (green), 5° (Hght blue), and 6° (blue) («y ^ 0). The 
higher frequencies become stronger with the increasing 
viewing angle. For Case A the power spectrum is scaled 
as P ^ 0)' as proposed for jitter radiation f§\. For all 
Cases the spectra are much steeper than the slope 1/3 
for the synchrotron radiation. 

The second row in Fig. 1 shows Case B with a larger 
jet velocity Vji 2 = 0.9924c with the other parameters 
kept the same as Case A. The spectra with larger view- 



ing angles are similar to those of Case A. The spectrum 
slope is smaller than that in Case A. However, due to the 
large jet velocity the higher frequencies at larger viewing 
angles (0°, 1°, 2°) become stronger On the other hand, 
with the smaller perpendicular velocities (Cases C and 
D), the gyroradius becomes very small. The spectra be- 
come weaker than those in Case A. As shown in the third 
and fourth rows in Fig. 1, the viewing angle dependence 
becomes very small. It should be noted that the slope 
of spectra is very steep for Case C. Spectral leakage is 
found. 

Cases D - F have a weaker magnetic fields (B^ = 
0.370) than Cases A - C. Case D has a small per- 
pendicular velocity. The trajectories are almost straight. 
The spectra look very similar to that for Bremsstrahlung 
[5]. The spectra become flat at lower frequencies. The 
peak spectral power is the weakest of all the cases. 
With larger perpendicular velocities the spectra become 
stronger than those with the smaller perpendicular veloc- 
ities. Case F shows the case with a larger time step (5 
times) with the same parameters as Case E. The spec- 
trum slope is very steep. The spectra in this case show 
two differences with those in Case E. First there exist 
positive slopes in the lower frequency. Second, due to 
the gyro-motion the spectra split due to the viewing an- 
gles. In particular the spectrum with larger viewing angle 
becomes stronger at high frequencies. 

As shown in Table 1, the critical angles for the off- 
axis radiation 0r = F^' are different. In this study we 
have obtained the off-axis radiation for the angles 0°, 
1°, 2°, 3°, 4°, 5°, and 6° (riy =^ 0). For cases D, E, and 
F the variation among different viewing angles is small 
since the angles are much smaller than the critical angle 
(25.3° and 13.35°). However, for case B (0r = 5°) the 



radiation shows larger differences for different viewing 
angles due to the small critical angle. For Case F (a 
longer time (340/cOpe) the spectra at high frequencies 
become stronger with larger viewing angles. 

These results validate the technique used in our code. 
It should be noted that the method based on the inte- 
gration of the retarded electric fields calculated by trac- 
ing many electrons described in the previous section can 
provide a proper spectrum in turbulent electromagnetic 
fields. On the other hand, if the formula for the frequency 
spectrum of radiation emitted by a relativistic charged 
particle in instantaneous circular motion is used [7., ,18], 
the complex particle accelerations and trajectories are 
not properly accounted for and the jitter radiation spec- 
trum is not properly obtained. 

DISCUSSION 

The procedure used to calculate jitter radiation using 
the technique described in the previous section has been 
implemented in our code. 

In order to obtain the spectrum of synchrotron (jit- 
ter) emission ll^SIlt], we consider an ensemble of elec- 
trons randomly selected in the region where the Ala- 
mentation (Weibel) instability |20] has fully developed, 
and electrons are accelerated in the generated magnetic 
fields. We calculate emission from about 20,000 elec- 
trons during the sampling time, fs = ^2 — h with Nyquist 
frequency cOn = 1 /2Ar where Af is the simulation time 
step and the frequency resolution Ao ~ 1/fs- However, 
since the emission coordinate frame for each particle is 
different, we accumulate radiation at fixed angles in sim- 
ulation system coordinates after transforming from the 
individual particle emission coordinate frame. This pro- 
vides an intensity spectrum as a function of angle rela- 
tive to the simulation frame x-axis (this can be any angle 
by changing the unit vector n in eq. (1)). A hypotheti- 
cal observer in the ambient medium (viewing the exter- 
nal GRB shock) views emission along the system jc-axis. 
This computation is carried out in the reference frame of 
the ambient medium in the numerical simulation. For an 
observer located outside the direction of bulk motion of 
the ambient medium, e.g., internal jet shocks in an ambi- 
ent medium moving with respect to the observer, an ad- 
ditional Lorentz transformation would be needed along 
the line of sight to the observer Spectra obtained from 
simulations can be rescaled to physical time scales. 

Emission obtained by the method described above is 
self-consistent, and automatically accounts for magnetic 
field structures on the small scales responsible for jitter 
emission. By performing such calculations for simula- 
tions with different parameters, we can then investigate 
and compare the quite contrasted regimes of jitter- and 
synchrotron-type emission (|9l[^[3|] for prompt and after- 
glow emission. The feasibility of this approach has been 



demonstrated and implemented 10, HI]. Thus, we will be 
able to address the issue of low frequency GRB spectral 
index violation of the synchrotron line of death [8]. 

Simulations incorporating jitter radiation are in 
progress using an MPI code |11] which speeds up 
considerably from the previous OpenMP code 1 17]. New 
results of jitter radiation will be presented separately. 
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